Characterization of age-associated gene expression changes in mouse sweat glands.

Evaporation of sweat on the skin surface is the major mechanism for dissipating heat in humans. The secretory capacity of sweat glands (SWGs) declines during aging, leading to heat intolerance in the elderly, but the mechanisms responsible for this decline are poorly understood. We investigated the molecular changes accompanying SWG aging in mice, where sweat tests confirmed a significant reduction of active SWGs in old mice relative to young mice. We first identified SWG-enriched mRNAs by comparing the skin transcriptome of Eda mutant Tabby male mice, which lack SWGs, with that of wild-type control mice by RNA-sequencing analysis. This comparison revealed 171 mRNAs enriched in SWGs, including 47 mRNAs encoding 'core secretory' proteins such as transcription factors, ion channels, ion transporters, and trans-synaptic signaling proteins. Among these, 28 SWG-enriched mRNAs showed significantly altered abundance in the aged male footpad skin, and 11 of them, including Foxa1, Best2, Chrm3, and Foxc1 mRNAs, were found in the 'core secretory' category. Consistent with the changes in mRNA expression levels, immunohistology revealed that higher numbers of secretory cells from old SWGs express the transcription factor FOXC1, the protein product of Foxc1 mRNA. In sum, our study identified mRNAs enriched in SWGs, including those that encode core secretory proteins, and altered abundance of these mRNAs and proteins with aging in mouse SWGs.

SUMOylation controls Hu antigen R posttranscriptional activity in liver cancer.

The posttranslational modification of proteins critically influences many biological processes and is a key mechanism that regulates the function of the RNA-binding protein Hu antigen R (HuR), a hub in liver cancer. Here, we show that HuR is SUMOylated in the tumor sections of patients with hepatocellular carcinoma in contrast to the surrounding tissue, as well as in human cell line and mouse models of the disease. SUMOylation of HuR promotes major cancer hallmarks, namely proliferation and invasion, whereas the absence of HuR SUMOylation results in a senescent phenotype with dysfunctional mitochondria and endoplasmic reticulum. Mechanistically, SUMOylation induces a structural rearrangement of the RNA recognition motifs that modulates HuR binding affinity to its target RNAs, further modifying the transcriptomic profile toward hepatic tumor progression. Overall, SUMOylation constitutes a mechanism of HuR regulation that could be potentially exploited as a therapeutic strategy for liver cancer.

UGDH promotes tumor-initiating cells and a fibroinflammatory tumor microenvironment in ovarian cancer.

BACKGROUND: Epithelial ovarian cancer (EOC) is a global health burden, with the poorest five-year survival rate of the gynecological malignancies due to diagnosis at advanced stage and high recurrence rate. Recurrence in EOC is driven by the survival of chemoresistant, stem-like tumor-initiating cells (TICs) that are supported by a complex extracellular matrix and immunosuppressive microenvironment. To target TICs to prevent recurrence, we identified genes critical for TIC viability from a whole genome siRNA screen. A top hit was the cancer-associated, proteoglycan subunit synthesis enzyme UDP-glucose dehydrogenase (UGDH). METHODS: Immunohistochemistry was used to characterize UGDH expression in histological and molecular subtypes of EOC. EOC cell lines were subtyped according to the molecular subtypes and the functional effects of modulating UGDH expression in vitro and in vivo in C1/Mesenchymal and C4/Differentiated subtype cell lines was examined. RESULTS: High UGDH expression was observed in high-grade serous ovarian cancers and a distinctive survival prognostic for UGDH expression was revealed when serous cancers were stratified by molecular subtype. High UGDH was associated with a poor prognosis in the C1/Mesenchymal subtype and low UGDH was associated with poor prognosis in the C4/Differentiated subtype. Knockdown of UGDH in the C1/mesenchymal molecular subtype reduced spheroid formation and viability and reduced the CD133 + /ALDH high TIC population. Conversely, overexpression of UGDH in the C4/Differentiated subtype reduced the TIC population. In co-culture models, UGDH expression in spheroids affected the gene expression of mesothelial cells causing changes to matrix remodeling proteins, and fibroblast collagen production. Inflammatory cytokine expression of spheroids was altered by UGDH expression. The effect of UGDH knockdown or overexpression in the C1/ Mesenchymal and C4/Differentiated subtypes respectively was tested on mouse intrabursal xenografts and showed dynamic changes to the tumor stroma. Knockdown of UGDH improved survival and reduced tumor burden in C1/Mesenchymal compared to controls. CONCLUSIONS: These data show that modulation of UGDH expression in ovarian cancer reveals distinct roles for UGDH in the C1/Mesenchymal and C4/Differentiated molecular subtypes of EOC, influencing the tumor microenvironmental composition. UGDH is a strong potential therapeutic target in TICs, for the treatment of EOC, particularly in patients with the mesenchymal molecular subtype.

Surfaceome analysis of extracellular vesicles from senescent cells uncovers uptake repressor DPP4.

Senescent cells are beneficial for repairing acute tissue damage, but they are harmful when they accumulate in tissues, as occurs with advancing age. Senescence-associated extracellular vesicles (S-EVs) can mediate cell-to-cell communication and export intracellular content to the microenvironment of aging tissues. Here, we studied the uptake of EVs from senescent cells (S-EVs) and proliferating cells (P-EVs) and found that P-EVs were readily taken up by proliferating cells (fibroblasts and cervical cancer cells) while S-EVs were not. We thus investigated the surface proteome (surfaceome) of P-EVs relative to S-EVs derived from cells that had reached senescence via replicative exhaustion, exposure to ionizing radiation, or treatment with etoposide. We found that relative to P-EVs, S-EVs from all senescence models were enriched in proteins DPP4, ANXA1, ANXA6, S10AB, AT1A1, and EPHB2. Among them, DPP4 was found to selectively prevent uptake by proliferating cells, as ectopic overexpression of DPP4 in HeLa cells rendered DPP4-expressing EVs that were no longer taken up by other proliferating cells. We propose that DPP4 on the surface of S-EVs makes these EVs refractory to internalization by proliferating cells, advancing our knowledge of the impact of senescent cells in aging-associated processes.

The YAP-TEAD complex promotes senescent cell survival by lowering endoplasmic reticulum stress.

Sublethal cell damage can trigger senescence, a complex adaptive program characterized by growth arrest, resistance to apoptosis and a senescence-associated secretory phenotype (SASP). Here, a whole-genome CRISPR knockout screen revealed that proteins in the YAP-TEAD pathway influenced senescent cell viability. Accordingly, treating senescent cells with a drug that inhibited this pathway, verteporfin (VPF), selectively triggered apoptotic cell death largely by derepressing DDIT4, which in turn inhibited mTOR. Reducing mTOR function in senescent cells diminished endoplasmic reticulum (ER) biogenesis, triggering ER stress and apoptosis due to high demands on ER function by the SASP. Importantly, VPF treatment decreased the numbers of senescent cells in the organs of old mice and mice exhibiting doxorubicin-induced senescence. Moreover, VPF treatment reduced immune cell infiltration and pro-fibrotic transforming growth factor-ß signaling in aging mouse lungs, improving tissue homeostasis. We present an alternative senolytic strategy that eliminates senescent cells by hindering ER activity required for SASP production.

Transcriptomes of human primary skin fibroblasts of healthy individuals reveal age-associated mRNAs and long noncoding RNAs.

Changes in the transcriptomes of human tissues with advancing age are poorly cataloged. Here, we sought to identify the coding and long noncoding RNAs present in cultured primary skin fibroblasts collected from 82 healthy individuals across a wide age spectrum (22-89 years old) who participated in the GESTALT (Genetic and Epigenetic Signatures of Translational Aging Laboratory Testing) study of the National Institute on Aging, NIH. Using high-throughput RNA sequencing and a linear regression model, we identified 1437 coding RNAs (mRNAs) and 1177 linear and circular long noncoding (lncRNAs) that were differentially abundant as a function of age. Gene set enrichment analysis (GSEA) revealed select transcription factors implicated in coordinating the transcription of subsets of differentially abundant mRNAs, while long noncoding RNA enrichment analysis (LncSEA) identified RNA-binding proteins predicted to participate in the age-associated lncRNA profiles. In summary, we report age-associated changes in the global transcriptome, coding and noncoding, from healthy human skin fibroblasts and propose that these transcripts may serve as biomarkers and therapeutic targets in aging skin.

NF-κB subunits direct kinetically distinct transcriptional cascades in antigen receptor-activated B cells.

The nuclear factor kappa B (NF-κB) family of transcription factors orchestrates signal-induced gene expression in diverse cell types. Cellular responses to NF-κB activation are regulated at the level of cell and signal specificity, as well as differential use of family members (subunit specificity). Here we used time-dependent multi-omics to investigate the selective functions of Rel and RelA, two closely related NF-κB proteins, in primary B lymphocytes activated via the B cell receptor. Despite large numbers of shared binding sites genome wide, Rel and RelA directed kinetically distinct cascades of gene expression in activated B cells. Single-cell RNA sequencing revealed marked heterogeneity of Rel- and RelA-specific responses, and sequential binding of these factors was not a major mechanism of protracted transcription. Moreover, nuclear co-expression of Rel and RelA led to functional antagonism between the factors. By rigorously identifying the target genes of each NF-κB subunit, these studies provide insights into exclusive functions of Rel and RelA in immunity and cancer.

DPP4 inhibition impairs senohemostasis to improve plaque stability in atherosclerotic mice.

Senescent vascular smooth muscle cells (VSMCs) accumulate in the vasculature with age and tissue damage and secrete factors that promote atherosclerotic plaque vulnerability and disease. Here, we report increased levels and activity of dipeptidyl peptidase 4 (DPP4), a serine protease, in senescent VSMCs. Analysis of the conditioned media from senescent VSMCs revealed a unique senescence-associated secretory phenotype (SASP) signature comprising many complement and coagulation factors; silencing or inhibiting DPP4 reduced these factors and increased cell death. Serum samples from persons with high risk for cardiovascular disease contained high levels of DPP4-regulated complement and coagulation factors. Importantly, DPP4 inhibition reduced senescent cell burden and coagulation and improved plaque stability, while single-cell resolution of senescent VSMCs reflected the senomorphic and senolytic effects of DPP4 inhibition in murine atherosclerosis. We propose that DPP4-regulated factors could be exploited therapeutically to reduce senescent cell function, reverse senohemostasis, and improve vascular disease.

Single-cell transcriptomic analysis uncovers diverse and dynamic senescent cell populations.

Senescence is a state of enduring growth arrest triggered by sublethal cell damage. Given that senescent cells actively secrete proinflammatory and matrix-remodeling proteins, their accumulation in tissues of older persons has been linked to many diseases of aging. Despite intense interest in identifying robust markers of senescence, the highly heterogeneous and dynamic nature of the senescent phenotype has made this task difficult. Here, we set out to comprehensively analyze the senescent transcriptome of human diploid fibroblasts at the individual-cell scale by performing single-cell RNA-sequencing analysis through two approaches. First, we characterized the different cell states in cultures undergoing senescence triggered by different stresses, and found distinct cell subpopulations that expressed mRNAs encoding proteins with roles in growth arrest, survival, and the secretory phenotype. Second, we characterized the dynamic changes in the transcriptomes of cells as they developed etoposide-induced senescence; by tracking cell transitions across this process, we found two different senescence programs that developed divergently, one in which cells expressed traditional senescence markers such as p16 (CDKN2A) mRNA, and another in which cells expressed long noncoding RNAs and splicing was dysregulated. Finally, we obtained evidence that the proliferation status at the time of senescence initiation affected the path of senescence, as determined based on the expressed RNAs. We propose that a deeper understanding of the transcriptomes during the progression of different senescent cell phenotypes will help develop more effective interventions directed at this detrimental cell population.

Pleiotropic effects of BAFF on the senescence-associated secretome and growth arrest.

Senescent cells release a variety of cytokines, proteases, and growth factors collectively known as the senescence-associated secretory phenotype (SASP). Sustained SASP contributes to a pattern of chronic inflammation associated with aging and implicated in many age-related diseases. Here, we investigated the expression and function of the immunomodulatory cytokine BAFF (B-cell activating factor; encoded by the TNFSF13B gene), a SASP protein, in multiple senescence models. We first characterized BAFF production across different senescence paradigms, including senescent human diploid fibroblasts (WI-38, IMR-90) and monocytic leukemia cells (THP-1), and tissues of mice induced to undergo senescence. We then identified IRF1 (interferon regulatory factor 1) as a transcription factor required for promoting TNFSF13B mRNA transcription in senescence. We discovered that suppressing BAFF production decreased the senescent phenotype of both fibroblasts and monocyte-like cells, reducing IL6 secretion and SA-ß-Gal staining. Importantly, however, the influence of BAFF on the senescence program was cell type-specific: in monocytes, BAFF promoted the early activation of NF-κB and general SASP secretion, while in fibroblasts, BAFF contributed to the production and function of TP53 (p53). We propose that BAFF is elevated across senescence models and is a potential target for senotherapy.

Transcription factor TFII-I fine tunes innate properties of B lymphocytes.

The ubiquitously expressed transcription factor TFII-I is a multifunctional protein with pleiotropic roles in gene regulation. TFII-I associated polymorphisms are implicated in Sjögren's syndrome and Lupus in humans and, germline deletion of the Gtf2i gene in mice leads to embryonic lethality. Here we report a unique role for TFII-I in homeostasis of innate properties of B lymphocytes. Loss of Gtf2i in murine B lineage cells leads to an alteration in transcriptome, chromatin landscape and associated transcription factor binding sites, which exhibits myeloid-like features and coincides with enhanced sensitivity to LPS induced gene expression. TFII-I deficient B cells also show increased switching to IgG3, a phenotype associated with inflammation. These results demonstrate a role for TFII-I in maintaining immune homeostasis and provide clues for GTF2I polymorphisms associated with B cell dominated autoimmune diseases in humans.

Enhanced myogenesis through lncFAM-mediated recruitment of HNRNPL to the MYBPC2 promoter.

The mammalian transcriptome comprises a vast family of long noncoding (lnc)RNAs implicated in physiologic processes such as myogenesis, through which muscle forms during embryonic development and regenerates in the adult. However, the specific molecular mechanisms by which lncRNAs regulate human myogenesis are poorly understood. Here, we identified a novel muscle-specific lncRNA, lncFAM71E1-2:2 (lncFAM), which increased robustly during early human myogenesis. Overexpression of lncFAM promoted differentiation of human myoblasts into myotubes, while silencing lncFAM suppressed this process. As lncFAM resides in the nucleus, chromatin isolation by RNA purification followed by mass spectrometry (ChIRP-MS) analysis was employed to identify the molecular mechanisms whereby it might promote myogenesis. Analysis of lncFAM-interacting proteins revealed that lncFAM recruited the RNA-binding protein HNRNPL to the promoter of MYBPC2, in turn increasing MYBPC2 mRNA transcription and enhancing production of the myogenic protein MYBPC2. These results highlight a mechanism whereby a novel ribonucleoprotein complex, lncFAM-HNRNPL, elevates MYBPC2 expression transcriptionally to promote myogenesis.

Transcriptomic analysis of human ALS skeletal muscle reveals a disease-specific pattern of dysregulated circRNAs.

Circular RNAs are abundant, covalently closed transcripts that arise in cells through back-splicing and display distinct expression patterns across cells and developmental stages. While their functions are largely unknown, their intrinsic stability has made them valuable biomarkers in many diseases. Here, we set out to examine circRNA patterns in amyotrophic lateral sclerosis (ALS). By RNA-sequencing analysis, we first identified circRNAs and linear RNAs that were differentially abundant in skeletal muscle biopsies from ALS compared to normal individuals. By RT-qPCR analysis, we confirmed that 8 circRNAs were significantly elevated and 10 were significantly reduced in ALS, while the linear mRNA counterparts, arising from shared precursor RNAs, generally did not change. Several of these circRNAs were also differentially abundant in motor neurons derived from human induced pluripotent stem cells (iPSCs) bearing ALS mutations, and across different disease stages in skeletal muscle from a mouse model of ALS (SOD1G93A). Interestingly, a subset of the circRNAs significantly elevated in ALS muscle biopsies were significantly reduced in the spinal cord samples from ALS patients and ALS (SOD1G93A) mice. In sum, we have identified differentially abundant circRNAs in ALS-relevant tissues (muscle and spinal cord) that could inform about neuromuscular molecular programs in ALS and guide the development of therapies.

Single-cell analysis of skeletal muscle macrophages reveals age-associated functional subpopulations.

Tissue-resident macrophages represent a group of highly responsive innate immune cells that acquire diverse functions by polarizing toward distinct subpopulations. The subpopulations of macrophages that reside in skeletal muscle (SKM) and their changes during aging are poorly characterized. By single-cell transcriptomic analysis with unsupervised clustering, we found 11 distinct macrophage clusters in male mouse SKM with enriched gene expression programs linked to reparative, proinflammatory, phagocytic, proliferative, and senescence-associated functions. Using a complementary classification, membrane markers LYVE1 and MHCII identified four macrophage subgroups: LYVE1-/MHCIIhi (M1-like, classically activated), LYVE1+/MHCIIlo (M2-like, alternatively activated), and two new subgroups, LYVE1+/MHCIIhi and LYVE1-/MHCIIlo. Notably, one new subgroup, LYVE1+/MHCIIhi, had traits of both M2 and M1 macrophages, while the other new subgroup, LYVE1-/MHCIIlo, displayed strong phagocytic capacity. Flow cytometric analysis validated the presence of the four macrophage subgroups in SKM and found that LYVE1- macrophages were more abundant than LYVE1+ macrophages in old SKM. A striking increase in proinflammatory markers (S100a8 and S100a9 mRNAs) and senescence-related markers (Gpnmb and Spp1 mRNAs) was evident in macrophage clusters from older mice. In sum, we have identified dynamically polarized SKM macrophages and propose that specific macrophage subpopulations contribute to the proinflammatory and senescent traits of old SKM.

A BDNF-TrkB autocrine loop enhances senescent cell viability.

Cellular senescence is characterized by cell cycle arrest, resistance to apoptosis, and a senescence-associated secretory phenotype (SASP) whereby cells secrete pro-inflammatory and tissue-remodeling factors. Given that the SASP exacerbates age-associated pathologies, some aging interventions aim at selectively eliminating senescent cells. In this study, a drug library screen uncovered TrkB (NTRK2) inhibitors capable of triggering apoptosis of several senescent, but not proliferating, human cells. Senescent cells expressed high levels of TrkB, which supported senescent cell viability, and secreted the TrkB ligand BDNF. The reduced viability of senescent cells after ablating BDNF signaling suggested an autocrine function for TrkB and BDNF, which activated ERK5 and elevated BCL2L2 levels, favoring senescent cell survival. Treatment with TrkB inhibitors reduced the accumulation of senescent cells in aged mouse organs. We propose that the activation of TrkB by SASP factor BDNF promotes cell survival and could be exploited therapeutically to reduce the senescent-cell burden.

Reduced endosomal microautophagy activity in aging associates with enhanced exocyst-mediated protein secretion.

Autophagy is essential for protein quality control and regulation of the functional proteome. Failure of autophagy pathways with age contributes to loss of proteostasis in aged organisms and accelerates the progression of age-related diseases. In this work, we show that activity of endosomal microautophagy (eMI), a selective type of autophagy occurring in late endosomes, declines with age and identify the sub-proteome affected by this loss of function. Proteomics of late endosomes from old mice revealed an aberrant glycation signature for Hsc70, the chaperone responsible for substrate targeting to eMI. Age-related Hsc70 glycation reduces its stability in late endosomes by favoring its organization into high molecular weight protein complexes and promoting its internalization/degradation inside late endosomes. Reduction of eMI with age associates with an increase in protein secretion, as late endosomes can release protein-loaded exosomes upon plasma membrane fusion. Our search for molecular mediators of the eMI/secretion switch identified the exocyst-RalA complex, known for its role in exocytosis, as a novel physiological eMI inhibitor that interacts with Hsc70 and acts directly at the late endosome membrane. This inhibitory function along with the higher exocyst-RalA complex levels detected in late endosomes from old mice could explain, at least in part, reduced eMI activity with age. Interaction of Hsc70 with components of the exocyst-RalA complex places this chaperone in the switch from eMI to secretion. Reduced intracellular degradation in favor of extracellular release of undegraded material with age may be relevant to the spreading of proteotoxicity associated with aging and progression of proteinopathies.

LncRNA OIP5-AS1-directed miR-7 degradation promotes MYMX production during human myogenesis.

Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) modulate gene expression programs in physiology and disease. Here, we report a noncoding RNA regulatory network that modulates myoblast fusion into multinucleated myotubes, a process that occurs during muscle development and muscle regeneration after injury. In early stages of human myogenesis, the levels of lncRNA OIP5-AS1 increased, while the levels of miR-7 decreased. Moreover, OIP5-AS1 bound and induced miR-7 decay via target RNA-directed miRNA decay; accordingly, loss of OIP5-AS1 attenuated, while antagonizing miR-7 accelerated, myotube formation. We found that the OIP5-AS1-mediated miR-7 degradation promoted myoblast fusion, as it derepressed the miR-7 target MYMX mRNA, which encodes the fusogenic protein myomixer (MYMX). Remarkably, an oligonucleotide site blocker interfered with the OIP5-AS1-directed miR-7 degradation, allowing miR-7 to accumulate, lowering MYMX production and suppressing myotube formation. These results highlight a mechanism whereby lncRNA OIP5-AS1-mediated miR-7 decay promotes myotube formation by stimulating a myogenic fusion program.

Single-Cell Analysis of the Transcriptome and Epigenome.

Epigenome regulation has emerged as an important mechanism for the maintenance of organ function in health and disease. Dissecting epigenomic alterations and resultant gene expression changes in single cells provides unprecedented resolution and insight into cellular diversity, modes of gene regulation, transcription factor dynamics and 3D genome organization. In this chapter, we summarize the transformative single-cell epigenomic technologies that have deepened our understanding of the fundamental principles of gene regulation. We provide a historical perspective of these methods, brief procedural outline with emphasis on the computational tools used to meaningfully dissect information. Our overall goal is to aid scientists using these technologies in their favorite system of interest.

Systematic Identification of circRNAs in Alzheimer's Disease.

Mammalian circRNAs are covalently closed circular RNAs often generated through backsplicing of precursor linear RNAs. Although their functions are largely unknown, they have been found to influence gene expression at different levels and in a wide range of biological processes. Here, we investigated if some circRNAs may be differentially abundant in Alzheimer's Disease (AD). We identified and analyzed publicly available RNA-sequencing data from the frontal lobe, temporal cortex, hippocampus, and plasma samples reported from persons with AD and persons who were cognitively normal, focusing on circRNAs shared across these datasets. We identified an overlap of significantly changed circRNAs among AD individuals in the various brain datasets, including circRNAs originating from genes strongly linked to AD pathology such as DOCK1, NTRK2, APC (implicated in synaptic plasticity and neuronal survival) and DGL1/SAP97, TRAPPC9, and KIF1B (implicated in vesicular traffic). We further predicted the presence of circRNA isoforms in AD using specialized statistical analysis packages to create approximations of entire circRNAs. We propose that the catalog of differentially abundant circRNAs can guide future investigation on the expression and splicing of the host transcripts, as well as the possible roles of these circRNAs in AD pathogenesis.